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Abstract 

For a given, arbitrary graph, what is the epidemic threshold? That is, under what conditions 
will a virus result in an epidemic? We provide the super-model theorem, which generalizes 
older results in two important, orthogonal dimensions. The theorem shows that (a) for a wide 
range of virus propagation models (VPM) that include all virus propagation models in standard 
literature (say, |8||5|), and (b) for any contact graph, the answer always depends on the first 
eigenvalue of the connectivity matrix. We give the proof of the theorem, arithmetic examples 
for popular VPMs, like flu (SIS), mumps (SIR), SIRS and more. We also show the implications 
of our discovery: easy (although sometimes counter-intuitive) answers to 'what-if questions; 
easier design and evaluation of immunization policies, and significantly faster agent-based 
simulations. 
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1 Introduction - Preliminaries 



Given a social or computer network, where the links represent who has the potential to infect 
whom, can we tell whether a virus will create an epidemic, as opposed to quickly becoming 
extinct? This is a fundamental question in epidemiology; intuitively, the answer should depend 
on (a) the graph connectivity and (b) the virus propagation model (VPM). This threshold is the 
minimum level of virulence to prevent a virus from dying out quickly [14 J. 

No result till now has unified the varied observations on the effect of different network struc- 
tures on the way different diseases spread. Moreover, with the exception of our earlier work on 
the SIS model [4] and its follow-up [6]; no other analysis examines arbitrary-topology graphs: 
the overwhelming majority of work focuses either on full-clique topologies (everybody contacts 
everybody else), or on 'homogeneous' graphs fl"2l [13] , or on power-law graphs [18] or hierar- 
chical (near-block-diagonal) topologies 1 9 1 (people within a community contact all others in this 
community, with a few cross-community contacts). 

We show that, irrespective of the virus propagation model, the effect of the underlying topol- 
ogy can be captured by just one parameter: the first eigenvalue Aj of the adjacency matrix A. 
In particular we cover all models given in the standard survey by Hethcote [8 J, which includes 
models like SIS (no immunity, like flu - 'susceptible, infected, susceptible') and SIR (life-time im- 
munity, like mumps: 'susceptible, infected, recovered'). We also include numerous other cases 
like SIRS [5 J (temporary immunity), our own useful generalizations SIV (vigilance/ vaccination 
with temporary immunity) and SEIV (vigilance / vaccination with temporary immunity and virus 
incubation) and many more. A few of these models are shown in Figure [l] organized in a lattice. 

Informally, our result can be stated as: 

Informal Theorem 1. For any virus propagation model (VPM) in the published literature, operating 
on an underlying undirected contact-network of any arbitrary topology with adjacency matrix A, the epi- 
demic threshold depends only on the first eigenvalue 

M 

of A and some constant Cvpm that is determined by the virus propagation model. 

Next, we give some preliminary definitions and in the upcoming sections, state our main re- 
sult together with its application on a select few standard models and its potential uses (Section|2] 
and Section |3| and finally give simulation results to illustrate the super-model theorem and dis- 
cuss some of its implications (Section |5|. The proof roadmap is explained in Section [4] with a 
description of the general model while the complete proof is presented in the Appendix. 

1.1 Preliminaries 

Table [T] and Table [2] list common terminology and describe some of the epidemic models we will 
be using in our article. Specific mathematical symbols and notations used for proofs are described 
in Appendix [A] 

Figure [T] shows the generalization hierarchy for some common epidemic models. The brown 
colored nodes denote standard VPMs found in literature while the blue colored nodes denote 
our generalizations. Each VPM is a generalization of all the models below it e.g. SIV is a gen- 
eralization of SIRS, SIR and SIS. Our main generalization, ST 2 V*, is described in detail later in 
Section|4j] 
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Figure 1: Virus Propagation model hierarchy (actually, lattice) for some standard models in- 
cluding SIRS (temporary immunity), SIV (vigilance, i.e., pro-active vaccination); SEIV (in- 
cludes the 'exposed but not infectious' state, and temporary vigilance); MSEIR (with the pas- 
sive immune state M); and our main generalization S*I 2 V*. The brown colored nodes denote 
standard VPMs found in literature while the blue colored nodes denote our generalizations. 
Each VPM is a generalization of all the models below it. 



3 



Table 1: Common Terminology 



VPM 


virus-propagation model 


NLDS 


non-linear discrete-time dynamical system 


P 


attack/ transmission probability over a contact-link 


6 


healing probability once infected 


y 


immunization-loss probability once recovered (in SIRS) or vigilant (in 
SIV, SEIV) 


£ 


virus-maturation probability once exposed - hence, 1 — e is the virus- 
incubation probability 


e 


direct-immunization probability when susceptible 


A 


adjacency matrix of the underlying undirected contact-network 


N 


number of nodes in the network 


Ai 


largest (in magnitude) eigenvalue of A 


s 


effective strength of a epidemic model on a graph with adjacency matrix 
A 


Table 2: Some Virus Propagation Models (VPMs) 


SIS 


'susceptible, infected, susceptible' VPM - no immunity, like flu 


SIR 


'susceptible, infected, recovered' VPM - life-time immunity, like mumps 


SIRS 


VPM with temporary immunity 


SIV 


'susceptible, infected, vigilant' VPM - immunization/ vigilance with tem- 
porary immunity 


SEIR 


'susceptible, exposed, infected, recovered' VPM - life-time immunity and 
virus incubation 


SEIV 


VPM with vigilance /immunization with temporary immunity and virus 
incubation 



2 Main Result 

The tipping point for each of the models captures a fundamental transition in the behavior of the 
system between the two extremes: a network-wide epidemic, versus a minor local disturbance 
that fizzles out quickly. We use the typical definition of the threshold used in literature [4.8.2. 6J. 
Intuitively below threshold the virus has no realistic chance of spreading the infection while 
above threshold the virus can take over and create an epidemic. For the SIS model, the tipping 
point describes some maximum strength of a virus, that will guarantee no epidemic [4, 8J. We 
define strength a little later (see Equation^. Similarly for the SIR model, the tipping point relates 
the explosiveness of the infection phase w.r.t. the virus strength, since in this model the virus will 
become extinct. 

In order to standardize the discussion of threshold results, we cast the threshold problem as 
expressing the normalized effective strength of a virus as a function of the particular propagation 
model and the particular underlying contact-network. So we are 'above threshold' when the 
effective strength s > 1, 'under threshold' when s < 1 and the threshold or the tipping point is 
reached when s = 1. 

Intuitively the effective strength s can be thought of as the basic reproduction number Rq 18 1 
frequently used in epidemiology. The effective strength is roughly the generalized Ro for the 
virus model and an arbitrary graph and is the quantity which determines the tipping point of an 
infection over a contact-network. 
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Formally, our main result is: 



Theorem 1 (Super-model theorem - sufficient condition for stability). For virus propagation models 



which satisfy our general initial assumptions (see Section 4.1 1 and for any arbitrary undirected graph with 



adjacency matrix A and largest eigenvalue \\, the sufficient condition for stability is given by: 

(1) 



s < 1 



where, s (the effective strength) is: 



s = Ai • Cvpm (2) 



and Cvpm is a constant dependent on the virus propagation model (given by Equation 30 in the Appendix). 
Hence, the tipping point is reached when s = 1. 

2.1 Two, Orthogonal Generalizations 

Our result generalizes along two different, difficult directions: (a) arbitrary contact-network 
topologies and (b) several virus propagation models (VPMs). 

2.1.1 General Topologies 

Much of previous work 1 1 ] has concentrated on the analysis of VPMs on specific types of contact- 
networks, typically cliques or homogeneous graphs. We include them all, as special cases. Specif- 
ically 

• Cliques, where every node contacts every other node. In that case, our result gives Aj = N, 
where N is the number of nodes in the graph 

• Homogenous graphs, with fixed degree d and random Erdos-Renyi graphs with expected 
degree d (e.g. see [12, 13|). In all these cases we have Ai = d, and our theorem includes the 
previous results. 

• Hierarchical (i.e., near-block-diagonal), e.g. |9j. 

• Power-law (e.g. I18|) 

Theorem [T] provides a simple and natural generalization of these results to arbitrary graphs. 
For example, previous results [18 1 have shown that the epidemic threshold in case of scale-free 
(power-law) networks is vanishingly small as the size N of the network increases. This is a 
corollary of our theorem: When a power-law graph grows (N — > oo), the largest eigenvalue 
grows with the highest degree, which also grows infinity, and thus the threshold approaches 
zero. 



2.1.2 General VPMs 



We also gen eralize with respect to various VPMs which satisfy our very general assumptions (see 



Section 



4.1 



i. We refer to our generalized model as ST 2 V, because it has an arbitrary number of 
susceptible states, two infectious / infected states, and an arbitrary number of vigilant / vaccinated 
(= recovered) states. All the standard models (like see [8], [5]) are simply special cases of S*I 2 V*: 

• the typical flu model, SIS, is a special case. 

• the typical mumps model, SIR, which corresponds to permanent immunity, is a special case, 
with one state for each class - S belongs to the Susceptible class, I belongs to the Infected 
class and R belongs to the Vigilant class. 

• the SIRS model (temporary immunity), similar to the SIR model 

• the SEIRS model ([8 ], page 601) where the virus has an incubation period (state 'E': exposed, 
but not infectious), and all other ingredients of the SIRS model (temporary immunity). 
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and more like SIV, SEIV, MSEIV etc. which generalize some specific models. The generalization 
hierarchy is shown in Figure]!] We elaborate more on the above in Section [472] 

We now give a brief summary of our threshold results (Table |3J by applying Theorem [T] on 
some standard epidemic models. Note the effect of the contact-network in effective strength for 
each model is captured solely by one parameter, ~K\ the first eigenvalue of the adjacency matrix 
of the network. Again, our result is a general one and these models just highlight the ready 
applicability of the result to standard VPMs in use. 



Table 3: Threshold results for some models. 



Model 


Effective Strength (s) 


Threshold (tipping point) 


SIS 


- (D 




SIR 


* (!) 




SIRS 


- (!) 




SIV 


Ar (5(7+6)) 


S = 1 


SEIR 


- (!) 




SEIV 


A i- (5(7+6)) 




SIiI 2 ViV 2 {used to model the H.I.V. 
virus, e.g. see 0) 


a ( Piv 2 +|3 2 eA 
Ar \v 2 (e+v 1 ) J 





3 Applications 

The results in this article can be fundamental to numerous applications. We describe a few im- 
portant ones next. 

3.1 Fast answers to "what-if " questions and guiding policy 

The threshold results can help quickly determine the result of plausible situations. For exam- 
ple, what happens if the virus is twice as infectious (virulent)? Similarly, what happens when 
there is a weaker strain of the virus? Our results will help in determining whether there is a 
danger of the infection taking-off or not. Naturally then this can feed into policy decisions for 
controlling epidemics. Assuming some models for the underlying contact network (like scale- 
free, small-world, hierarchical etc.) we can estimate which nodes/classes should be quarantined 
or immunized first. Given the linear dependence on Ai, we want to immunize nodes (and hence 
remove them from the contact graph) which will drop the Aj value the most so that the resultant 
infection becomes below threshold and dies out. For example, they may decide to immunize 
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teachers and kindergarten children first to control the epidemic. In addition, they can impose re- 
strictions on travel so as to not increase the Aj and hence the effective strength for the virus. The 
above discussion also illustrates the generality of our result. Policy makers can assume any graph 
model which captures the contact behavior of the population the best and still use our threshold 
result to guide policy. 

A lot of work has been done to show that immunizing high-degree nodes in scale-free net- 
works is a good idea because of the vanishing threshold result IfTHl . But significantly just con- 
centrating on high-degree nodes will miss those low-degree nodes which are good "bridges" 
and hence can have a important influence on decreasing Ai when immunized. Intuitively how 
disparate the groups are to which a node connects is also important in addition to how many 
groups one is connected to. For example, a single common friend of only some sportspeople and 
movie stars can have a huge impact in the outbreak of a disease even if he / she knows only a few 
sportspeople and movie stars (while sportspeople and movie stars are themselves very tightly 
connected). 

We have been concentrating on biological viruses only. But various biological virus models 
have been used to model computer viruses as well [14] e.g. J?] introduced the SHIR model ('sus- 
ceptible', 'hidden', 'infected', 'recoverable') to model computers under email attack. More so 
than the biological cases, it is easier to get the entire underlying network. Hence our threshold 
results can be precisely used to make the network more robust to malware and computer viruses 
by selectively "removing" nodes from the contact-network by immunizing them like installing a 
firewall on them etc. 

3.2 Simplifying epidemiological simulations 

Epidemiological simulations in general have several parameters and are computationally expen- 
sive to run. By distilling the impact of the network topology and also giving the exact threshold, 
we can greatly simplify these efforts. For example, parameters which do not affect the effective 
strength of the contagion need not be varied. In addition, regular topologies like cliques, block- 
diagonal matrices lend themselves to fast eigenvalue computations. This can be taken advantage 
of to quickly identify parameter spaces where simulations would be useful. 

3.3 Viral Marketing 

A variety of dynamic processes on graphs are modeled like epidemic spreading. In contrast to 
the biological viruses, conversely, we may actually want the spread of a contagion as quickly as 
possible in some situations e.g. spread of a product or idea in a network of individuals. The 
Bass model [3] fits product adoption data using parameters for pricing and marketing effects 
ignoring topology and hence assuming that all adopters have equal probability of influencing 
non-adopters. A more refined picture using our result can be constructed of when a product 
gains massive adoption on a network (equivalent to an "epidemic"). Marketers can also then 
decide which target groups should the product be more marketed to and how aggressive should 
the marketing be for an estimated network of individuals. 

3.4 Software patch distribution and more 

Another application is efficient spreading of software patches over a computer network. The 
patches behave like computer worms |19| and can help defend against other malicious worms. 
We want to maximize the spreading of the patch over the network. Given full knowledge of the 
router-network involved, we can estimate how "aggressive" the patch-worm (say by increasing 
the number of probes for possible hosts before dying out) has to be to initiate an "epidemic". Our 
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threshold results can also help determine the "vulnerability" and hence consequently the cost of 
not patching a part of the network say after one wave to "worm-patching". 

Various epidemic models have also been used to model blog cascades. These models can be 
now applied to arbitrary graphs e.g. study propagation of memes through blogs 1 15 1. 

4 Proof Roadmap 

The first basic idea behind the proof is approximating the epidemiological models by a discrete 
time non-linear dynamical system (NLDS). A NLDS can be represented by 

Pt+i = g(P t ) 

where g is some non-linear function operating on a vector. We define the vector P t such that it 
specifies the state of the system at time t; the exact definition will differ from model to model 
but it effectively encodes the probability of each node in the graph of being in any given state at 
time t. An equilibrium point (also called a fixed point) of the NLDS is the state vector (i.e. some 
particular P) of the NLDS which does not change. Thus at the equilibrium point P t +i = Pt = x. 
Formulating a NLDS helps us in leveraging the vast results in NLDS literature. In particular, 
Theorem [2] (given in Appendix |Bj gives the conditions when an equilibrium point a NLDS is 
stable (Figure [2] illustrates this concept). It relates the eigenvalues of the Jacobian of the NLDS at 
the equilibrium point with the stability of the NLDS at that point. 




(A) Unstable (B) Stable (C) Neutral (at threshold) 

Figure 2: Different types of equilibria for a discrete-time NLDS: Unstable, Stable and Neutral. 
The equilibrium is neutral at threshold. 

Intuitively, the tipping point for any model then deals with analyzing the stability of the 
corresponding NLDS at the point when none of the nodes in the graph are infected, because 
otherwise the infection can still spread. If the equilibrium is unstable, a small "perturbation" 
(physically in the form of a few initial nodes getting infected) will push the system further away 
which physically means more and more nodes will get infected leading to an epidemic. But if the 
equilibrium is stable, the system will try to come back to the fixed point without going "too-far" 
away in effect, "controlling the damage". At threshold, the tendencies to go further away and 
come-back will be the same. In other words, the equilibrium is stable below the threshold and is 
neutral at the tipping point. 

The requirement imposed by Theorem [2] on the eigenvalues of the 'Jacobian' of the corre- 
sponding NLDS for any virus propagation model eventually reduces to a simple condition on 
the eigenvalue of the adjacency matrix. This condition translates into the effective strength of the 
virus under the model. The reason we can reduce the condition to one on the adjacency matrix is 
due to the special structure of the models, which is explained next. 

4.1 General Model 

Suppose the virus-propagation model has m states (e.g. m = 3 for the SIR model with states S, I 
and R) and it operates on a graph of N nodes. Consider then a m. • N x 1 vector P t which captures 
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Endogenous 
Transitions 




Figure 3: General State Diagram for a node in the graph - it is not a simple Markov chain. In 
essence, there are only three types of state 'classes' - Susceptible, Infected and Vigilant. Only 
cross-class transitions have been shown. Note the unidirectional arrow from Infected to Vigi- 
lant. Only one type of transition is exogenous i.e. graph-based (red curvy arrow) affected only 
by the neighbors of the node, all other transitions shown are endogenous. Any transitions 
between the states of the same class are also endogenous. The red arrow always ends at I lr 
i.e. any state in the Infected class can cause a graph-based transition from any state in the 
Susceptible class only to the Ii state. 



the probability of each node being in any of m states at a given time t. Specifically: 



P Sl/ l,t 

Psi,2,t 

P[ - P Sl ,N,t (3) 
Ps 2 ,l,t 

s m ,N,t 

where, P s . xt — the probability of node i being in state s j at time t. This vector completely defines 
the system at time t and our NLDS equation will track the evolution of this vector across time. 

Every epidemic model will have some fundamental states and the choice of which states to 
include in a model depends on the particular disease characteristics |8J. In fact, every model can 
be essentially thought of having states in any of the following broad classes (see Figure [3b: 

Susceptible Class: Nodes in any of the states in this class signify that it (the individual) is sus- 
ceptible and can get infected by any neighboring node in a state of the Infected class. 

Infected Class: In a state of this class, the node is infectious in the sense that it is capable of trans- 
mitting the infection to its neighbors. Note that each such state will have a transmissibility 
parameter (e.g. (3 in the SIR model for the infectious state I). This definition is general 
enough to accommodate models with states which have their transmissibility parameter = 
i.e. they are 'exposed' but not infectious (e.g. the E state in the SEIR model is a state which 
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is in the Infected class in the sense that it can potentially cause infections but is not by itself 
infectious). 

Vigilant/Vaccinated Class: This class contains states which are not in either of the other two 
classes Susceptible and Infected. Nodes in any of the states in this class can not get infected 
nor do can they potentially cause infections. States like M (called as the passive immune 
state which feeds to the susceptible state), R (the recovered/died state where the node either 
gets permanent immunity or dies and hence doesn't participate in the epidemic further) etc. 
are all conceptually of the Vigilant type. Notice that importantly by definition, any state in 
the Vigilant class will not have a direct transition to any state in the Infected class. This is so 
because otherwise that state can be potentially infectious and thus is not part of the Vigilant 
class. 

Apart from the types of states, models also fundamentally have only two types of transitions: 
Exogenous (graph-based, in particular affected only by the neighbors) and Endogenous (caused by 
the node itself by some probability at every time step). For example, the transition from S to I 
in the SIRS model is an exogenous transition while the transition from R to S is an endogenous 
transition. It is the presence of the graph-based transitions that makes our model not a simple 
Markov chain and brings in the topology of the graph into play. In our general model note that 
only one class of transitions is graph-based, all others are endogenous transitions. The exogenous 
transitions are the ones that take a node from a state in the Susceptible class to the Infected class. 
In addition we assume any such exogenous transition ends at the Ij state; but these transitions 
will themselves be caused by some neighbor node in any of the states of the Infected class. 

Assumption 1 (Transition Assumption 1). The only way to get infected is through your neighbors i.e. 
there is no path to a state in the Infected class from a state in the Susceptible class composed solely of 
endogenous transitions. 

Assumption 2 (Transition Assumption 2). Any exogenous (graph-based) transition always results in 
a transition from a state in the Susceptible class to the l\ state. 

As seen from a single node, these models will look like in Figure [3] It gives the general 
state diagram for a node in the graph together with the assumptions on the transitions discussed 
above. We have shown only cross-class transitions and their types. All transitions between states 
of the same class are of course endogenous. 

4.2 Examples 

Our general characterization is powerful enough to seamlessly capture all the practical mod- 
els like SIS, SIR, SIRS, SIV, SEIR, SERIS, MSIR, MSEIR etc. EJ while still being mathematically 
tractable to yield simple threshold equations. 




Figure 4: Transition diagram for the SEIV model. 
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We give one instantiation of our general model in the form of the SEIV model, which itself 
generalizes several known VPMs (see Figure III. SEIV has one state (S) in the Susceptible class, 
two (E and I) in the Infected class and one (V) state in the Vigilant class. The transition diagram 
for the SEIV model in shown in Figure [4] Note the similarities in Figures [3] and [4] e.g. E state is 
the Ii state of Figure [3] The infection is caused only by the I state and the E is only a latent state. 
The SEIR, SEIRS, SIRS, SIV, SIR and SIS models are all special cases of SEIV: 

• SIS, is a special case, with e = l,y = 1, = 0. 

• SIR, which corresponds to permanent immunity, with e — l,y — 0, 9=0. 

• SIRS (temporary immunity), similar to the SIR model, with e = 1, 8 = 0. 

• SIV, our own model, where ' V stands for vigilant or vaccinated, with e = 1 . 

• SEIR, similar to SIR but where the virus has an incubation period, with y = 0, 8 =0. 

• SEIRS where the virus has an incubation period, and all other ingredients of the SIRS model, 
with 8=0. 

The complete proof for the general model with two states in the infected class (we call this as the 
S*I 2 V*generalized model) is given in the Appendix. The proof follows the following roadmap: 

1. We develop the system equations of our general NLDS in Appendix [B] 

2. The relevant fixed point is computed in Appendix |C| 

3. The Jacobian at the fixed point is constructed in Appendix |D| 

4. Finally we prove Theorem[T]by finding the condition to bound the eigenvalues of the Jaco- 
bian in Appendix|E]and Appendix|F] 

5 Discussion 

We discuss some simulation examples for our threshold in some models and a few direct impli- 
cations of the super-model theorem in this section. We also illustrate what the result implies for 
the 'vulnerability' of the underlying contact graph for epidemics. Apart from the dependence of 
the threshold on Ai, it is instructive to note some unexpected results in specific models as well. 

5.1 Simulation Examples 

We conducted some computer simulation experiments on the Oregon AS router graprj^to illus- 
trate our super-model theorem. Figures [5j [6] [7] and [8] give an overview of the simulations for the 
SIR, SEIR, SIRS and SEIV models. All values are average over several runs of the simulations. 
In short, as expected from the theorem, the difference in behavior above, below and at threshold 
can be distinctly seen for each of these models. 

Figures |5]A),[6]A),[7]A) and|8jA) show a time-evolution plot of the fraction of infected nodes 
in the graph for different values of the effective strength of the virus specifically above threshold (in 
red) and under threshold (in green). Note the qualitative difference in the behavior of the system 
above and below the threshold. 

Figure |5| A) and Figure |6jA) deal with SIR and SEIR which don't have a steady state because 
of the recovered state. The number of susceptibles available in the graph decrease with each new 
infection and hence epidemics will disappear unlike in models like SIS and SIRS. Additionally 
as a result, above the threshold, the number of infections has an explosive phase and then they 
go down to zero while in the SIRS and SEIV models (see Figures [7tA) and|8|A)) the number of 
infections reach a steady state value. Contrast this with the under threshold behavior, where in 
all the models the number of infections aggresively go down to zero. 

1 This is a real network graph collected from the Oregon router views. It contains 15,420 links among 3,995 AS peers. 
More information can be found fromlhttp : / /topology . eecs .umich . edu/data . htmll 
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(A) SIR Infective Fraction Time Plot (log-log) 



| 1000 



Our Threshold 
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Effective Strength of Virus 

(B) SIR Footprint vs Strength (lin-log) 



Figure 5: SIR (all values averages over several runs): (A) Plot of Infective Fraction of Pop- 
ulation vs Time (log-log). Note the qualitative difference in behavior under and above the 
threshold. (B) Plot of Final number of cured nodes (the footprint) vs Effective Strength (lin- 
log). Note the tipping point is exactly when the effective strength s = 1. 




Figure 6: SEIR (all values averages over several runs): (A) Plot of Infective Fraction of Popula- 
tion vs Time (log-log). Note the qualitative difference in behavior under and above the thresh- 
old. Also notice the initial "silent" period for above threshold because of virus-incubation. 
(B) Plot of Final number of cured nodes (the footprint) vs Effective Strength (lin-log). Note 
the tipping point is exactly when the effective strength s = 1. 
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SIRS Threshold 




Time ticks 

(A) SIRS Infective Fraction Time Plot (log-log) 



Our Threshold 
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(B) SIRS Max. Infections till steady 
state vs Strength (lin-log) 



Figure 7: SIRS (all values averages over several runs): (A) Plot of Infective Fraction of Pop- 
ulation vs Time (log-log). Note the qualitative difference in behavior under and above the 
threshold. (B) Plot of Max. number of infected nodes till steady state vs Effective Strength 
(lin-log). Note the tipping point is exactly when the effective strength s = 1. 
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(A) SEIV Infective Fraction Time Plot (log-log) 
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(B) SEIV Max. Infections till steady 
state vs Strength (lin-log) 



Figure 8: SEIV (all values averages over several runs): (A) Plot of Infective Fraction of Popula- 
tion vs Time (log-log). Note the qualitative difference in behavior under and above the thresh- 
old. Also notice the initial "silent" period for above threshold because of virus-incubation. 
(B) Plot of Max. number of infected nodes till steady state vs Effective Strength (lin-log). Note 
the tipping point is exactly when the effective strength s = 1. 
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Also note the initial "flat" period in the time plots for above threshold for the models having 
the Exposed (E) state, SEIR and SEIV. This is due to the virus-incubation period because of which 
there is an initial delay in number of infected nodes. This then results in an initial "silent" period 
after which the epidemic takes-off . As there is no such incubation period in SIR and SIRS, their 
plots do not show such silent periods. 

We also give a "take-off" plot for each model ((Figures gB), ^B), [TfB) and^B)). They show 
the final number of cured nodes (for SIR and SEIR) and max. number of infections till steady 
state (for SIRS and SEIV) vs the different strengths of the virus. Intuitively, these metrics measure 
the "footprint" of each infection. If the infection resulted in an epidemic then the footprint will 
be large. As our theorem predicted, the plots clearly illustrate that the tipping point in all the 
cases is at the point when the effective strength s = 1. 



5.2 Ai: Measure of Connectivity 



ooooo 

(A) Chain (M = 1.73) (B) Star (M = 2) (C) Clique (M = 4) 

Figure 9: Changing connectivity and vulnerability of graphs with changing Ai. Our result says 
that a virus has the most effective strength in a clique, star and chain, in that order. Thus the 
Clique is the most vulnerable. 

What does exactly the result mean w.r.t. the graph? Intuitively, Aj (also known as the spectral 
radius) of a graph captures the connectivity of the graph. The preceding threshold results show 
that this is precisely what matters in the epidemic threshold. Hence, more connected the graph 
is, more vulnerable it is to an epidemic by a virus. For example, see Figure [9] It shows three 
graphs with the same number of nodes in increasing Ai value. Note that although the star and 
chain (Figure [9] (A) and (B)) have the same number of edges, the star has a higher Aj value, thus 
boosting the effective strength s of a virus making the graph more vulnerable to epidemics. Also, 
recall that a d-regular graph (a homogeneous graph with all nodes of same degree d) has Ai = d. 
Hence, we can make the following observation: 

Observation 1 (Impact of Ai). Our threshold results suggest that an arbitrary graph behaves in the same 
way to a \\-regular graph. In other words, Ai captures the average neighborhood for a node in the graph. 

Notice that the actual dynamics of the epidemic may not be captured by Ai completely, but 
the threshold is solely dependent on the Ai (apart from the VPM). 




5.3 Insensitivity to Virus Incubation 

For the SEIV model which is similar to SIV, except for the E state, we have the following: 

Observation 2 (SEIV). The threshold here does not depend on the virus-maturation probability e i.e. the 
probability of transition from Exposed to Infected. 

This implies that the incubation period of the virus does not have an effect for purposes of 
the epidemic threshold. The parameter e, in effect, only delays / speeds-up the achievement of 



the threshold, not what the threshold itself is. As an example, this is also illustrated in Figure 10 
for the standard SEIR model |8| (it is also one of the special cases of our SEIV generalization). 
It shows the final number of cured nodes (the 'footprint') of computer simulations of the entire 
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SEIR Threshold - Effect of e 
2500 1 ■ — ■ — ■ — 




-2-1012 
10 10 1Q 10 10 

Effective Strength of Virus (X p/S) 

Figure 10: Changing e in the SEIR model on the AS graph (lin-log scale). Note that it does 
not affect the threshold (the tipping point is still at effective strength ^ = 1). All values are 
averages over multiple runs. 



virus infection as a function of the effective strength of the virus (Ai(3/6) on the AS graph. Note 
that the tipping point is the same for all the three colored curves which correspond to different 
values of e. 

5.4 Confirming "Prevention is better than cure" 

Clearly SIV is the most general model assuming no intermediate Exposed state and comprising 
three states (see Figure [Tj. We know its effective strength (Table [3} is: 

S = Ai • 



5(y + 6) 



By setting y, 8 etc. to zero we can get different models such as SIRS, SIR, SIS and so on. For 
example, when we set 9=0, the model reduces to SIRS, where the threshold we know is Aj (3/6. 
This generalization gives us further insight: 

Observation 3 (Rate of loss of immunity). Lowering the rate of loss of immunity i.e. having a smaller 
y (say due to better hygiene) decreases the effective strength s (and makes it harder for the virus to cause 
an epidemic) only so long as there is a mechanism to give a node direct immunity i.e. having a non-zero 
8 (say by using a vaccine) before an infection (in the Susceptible state) instead of after (in the Recovered 
state). 

Satisfyingly this fits well with the old adage 'Prevention is better than Cure'. This can also be 
vividly seen specifically in the case of SIRS (8 = 0): 

Observation 4 (SIRS). The threshold in SIRS does not depend on y, the immunity-loss rate, that is, the 
probability of becoming Susceptible again from the Recovered state. 



Figure 1 1 demonstrates this on the AS graph. Again, the tipping point does not depend on y, 
it is the same for all the three colored curves. This means that the rate of lapse or loss of immunity 
to a virus is immaterial to the threshold here. 
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SIRS Threshold - Effect of y 




-2-10 1 



'° Effective Strength of Virus (X 10 

Figure 11: Changing y in the SIRS model on the AS graph (lin-log scale). Note that it does 
not affect the threshold (the tipping point is still at effective strength ^ = 1). All values are 
averages over multiple runs. 

6 Conclusion 

In this article we provided two, orthogonal generalizations of earlier epidemic threshold results. 

• In the first direction, our result gives the threshold of the generalized S*I 2 V* model, which 
encompasses any epidemic model in published literature (|8| |5] etc.). 

• In the second direction, topology, we showed that for any, arbitrary, undirected contact- 
network, the effect of the topology can be captured solely by \, the first eigenvalue of the 
adjacency matrix. 

We also demonstrated simulation results which illustrate our result. In addition, we discussed 
some important applications and implications of our result for policy makers, scientists etc. like: 

• Fast answers to "what-if " questions. 

• Guiding immunization policies. 

• Simplifying epidemiological simulations. 

• Modeling Viral Marketing, Software patch distribution, Blog propagation etc. 
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Appendix 



A Notation 

Recall that we are dealing with the S*I 2 V* generalized model - it has two states Ii and I2 in the 
Infected class. To simplify notation, we refer to state Ii as E (the 'infection entrance state') and 
I2 as the I state in the proofs. The state E has a transmission probability of |3j and the state I has 
a transmission probability of p2- The states E and I here should be thought as to mean general 
infected states of our model and not in the sense of the specific E and I states in epidemic models 
like SEIR, SEIV etc. We also refer to the exogenous transitions as graph-based and endogenous 
transitions as internal interchangeably. Table [4] gives some of the notation we will be using in our 
description of the proof. 



Table 4: Notation and Symbols used in proofs 



m 


total number of states in the model 


q 


total number of states in the Susceptible and Vigilant classes of the model; 
hence m = q + 2 


w 


total number of states in the Susceptible class of the model 


Si, S2, . . . , s w 


general states in the Susceptible class 


E,I 


general states in the Infected class 


<Xku 


probability (constant and given) of transition from state K to state U 


Pi 


transmission probability for state E 




transmission probability for state I 


iu(E,I) 


probability that a node i does not receive any infections from E and I at 
time t 


X 


the fixed point vector our NLDS corresponding to when no node is in any 
of the Infected class states 




(same for each node) probability of being present in the S y state at x 


a 


lacobian matrix of the NLDS computed at x 



B System Equations 

We can develop the system equations i.e. explicitly specify the non-linear function g for the 



we assume 



NLDS based on the transition diagram of the model. As stated earlier in Section 4.1 
that infections are received only from infected neighbors i.e. those in states E and I the Infected 
class of states. Firstly, let's calculate the probability that a node i does not receive any infections 
in the next time step (call it Ci,t(E, I), E, I denotes that an infection is passed only from a neighbor 
in the E or I states). No infections are transmitted if: 

• Either a neighbor is not any of the infected states E and I 

• Or it is in state E and the transmission fails with probability 1 — (3i 

• Or it is in state I and the transmission fails with probability 1 — [3 2 

Since we assume infinitesimally small time steps (At — > 0), multiple events can be ignored for 
first-order effects in the time step. Also, assuming the neighbors are independent, we get: 

C U (E,I) = Yl (PE,j,t(l-|3i) + Pi,j,t(l-p 2 ) + (l-PE,i,t-Pi,j,t)) 

= Yl (l-A i , J (PiP E , J> t + P2Py,t)) ( 4 ) 
je{i..N} 
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where N£ (i) is the set of neighbors of node i in the graph. 

Also, the sum of probabilities of being in all the possible states for each node i should equal 
1. Hence, 

V u X Pk <u = 1 (5) 

K 

We can now write down the system equations as follows. A node i will be in any particular 
state S y of the Susceptible class at time t + 1 if : 

• Either it was in S y at time t and stayed in state S y i.e. it did not receive any infections from 
its neighbors and it did not change state internally from S y to any other state 

• Or it was in some other state U and changed state internally from U to Sy 

Hence, the probability of node i being in S y where S y is any state in the Susceptible class at 
time t + 1 is: 

Vi)=l,2 ( ...,w P Su ,i,t+i = X «KS s Put + Ps s ,u Ki,t(E,I)- Y- as « K ( 6 ) 
Similarly for the E state: 

w 

Pe,u+i= Y- a KE PK,i,t + X p s„u(l-Cu(E,I)) (7) 

K^S^Sz,...^ y=l 

and for any other state U ^ {Si, S2, . . . , S w , E}: 

Pu,i,t+i = Y k kuPk,u (8) 



As discussed earlier (Equation|3J, we can now define a probability vector P t by "stacking" all 
these probabilities which will completely describe the system at any time t and evolve according 
to the above equations. Note that the above equations are non-linear and naturally define the 
function g for the NLDS P t+1 = g(P t ). 

We have the following theorem about NLDS stability at a fixed point: 

Theorem 2 (Asymptotic Stability, e.g. see llTOl ). The system given by P t+1 = g(P t ) is asymptotically 
stable at an equilibrium point P = x, if the eigenvalues of3 — V9(*) are ' ess ^ nan 2 m absolute value, 
where, 

hi = [vg(x)] y = ^l?=x 

Hence, next we compute the fixed point we are interested in and the Jacobian of our NLDS at 
that point. 



C Fixed point 

We are interested in the stability of the equilibrium point (i.e. where Pt+i = Pt (= x)) of the NLDS 
which corresponds to when no one is infected. Only the transition from the Susceptible class 
states towards the Infected class states are graph-based (and can happen only when at least one 
of the nodes is in any of the Infected states), so the state-diagram for each node will be a simple 
Markov chain (call it MCsv) consisting of the Susceptible and Vigilant states (see Figure [12} . 
Note now there are no graph-based effects, hence each node is independent of others and will 
converge to steady state probabilities corresponding to the Markov chain. The steady state vector 
7t* (size q x 1, where q is the number of states in the Susceptible and Vigilant classes) which will 
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Figure 12: State Diagram for any node in the graph at the fixed point when no node is present 
in a state in the Infected class. Only cross-class edges are shown. Note that it is now a simple 
Markov chain with a unique steady state probability. 



be the same for each node can be computed from the following equations from standard Markov 
chain analysis: 



7t Tran 



MC S 



= 7T* & 



(9) 



i=l 



Hence n* is a probability vector and is the left eigenvector corresponding to eigenvalue 1 of 
the stochastic matrix TranMc sv or tne Markov chain MCsv- The full (m x 1) probability vector 
p~* for each node at this steady state will have the entries in n* for states in the Susceptible and 
Vigilant classes and for all states in the Infected class. The fixed point of the global original 
NLDS x can be finally represented as: 

[V* 

(10) 



where p~* is repeated N times (once for each node in the graph). Let pg be the steady state 
probability value in the vector p~* corresponding to the S y state. In other words, each node will 
have a probability of p£ of being present in the S y state at the fixed point. Also define, 



Ps 



(u) 



i.e. Pg is the total probability of each node at the fixed point of being present in any of the states 
of the Susceptible class. 



D The Jacobian 

We know from Theorem [2] that x is stable if the eigenvalues of 3 — V9(*0 are ^ ess than 1 in 
absolute value. From the definition of 3 we can see that it is a m N x m • N matrix with m (for 
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each state) square blocks of size N x N each (corresponding to every node in the graph). We can 
calculate 2 to be (states have been mentioned on the top and side for ease of exposition and I is 
the identity matrix of size N x N): 





s y 


K 




E 


I 




P-Ik^.e «s 1( k)I 






















u 


as y ul 


cxkuI 




a EU I 


ami 














E 




cxkeI 




a EE I + PsPiA 


ai E I + Psp2A 


I 


as, i 1 


a K iI 






anl 



where K is any state 7^ {E, 1} and U is any state 7^ {Si, S2, • • ■ , S w , E, I}. 

Recall the properties we are assuming for the epidemic models discussed in Section |4~T| 
see Figure[3j. Crucially they imply Vk=^e,i «-ke = and Vk^e,i <xki = 0- Hence 8 reduces to: 





s y 


K 




E 


I 










aES.I-Ps.PiA 


0C is li I-PsP2A 














u 




a KU I 






ot m I 














E 


0n,n 


0n,n 




cxeeI + pIPiA 


aiEl + p5P 2 A 


I 


0n,n 


0n,n 




cxeiI 


ani 



where On n is a N x N matrix with all zeros. 



E Eigenvalues of the Jacobian 



Note that d is very structured and can be written as: 

Bi 

02N,(m-2)N 



B 2 
B 3 



(12) 



where Bj, B 2 and B 3 are matrices of size (m - 2)N X (m - 2)N, (m - 2)N x 2N and 2N x 2N 
respectively. B3 corresponds to the E and I rows and columns of d i.e.: 



<xeeI + PsI3iA oci E I + Psp2A' 



(13) 



Bi and B 2 are defined similarly. Consider any eigenvector v (size mN x 1) and corresponding 
eigenvalue Ag of 3- We can write v as being composed of vector vi of size (m — 2)N x 1 and vector 
v 2 of size 2N x 1 i.e: 



v 2 



(14) 
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Also x and Ag satisfy the eigenvalue equation: 

0V = AgV 

: we get: 



Substituting from Equations 12 and 
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B 3 



Equation [16] implies the following the two relations: 



B1V1 + B 2 v 2 
B 3 v 2 



AgVi 
Agv 2 



(15) 



(16) 



(17) 
(18) 



From Equation 18 we can infer that precisely one of the following holds: 



1. v 2 =0 

2. v 2 is the eigenvector of B 3 (and consequently Ag is the matching eigenvalue of B 3 ) 



If v 2 = 0, Equation 17 reduces to 



BjVi = AgVi 

wherein again, either vi = or Ag is an eigenvalue of Bi. The condition vi = 6 is not meaningful 
as then v = (v is an eigenvector of 3 implies v is non-zero). Therefore the eigenvalues of 3 are 
given by the eigenvalues of Bi (with v 2 = 6) and the eigenvalues of B3. 



E.l Eigenvalues of 

From the expression for 3 derived in Section[D] note that: 

B x = T<g>I 

where <S> is the Kronecker product of two matrices and 

T 



(19) 



^s„u 



(20) 



We know from matrix algebra [11] that if C = E then Ca = D A ®E A , where C\ denotes a diag- 
onal matrix with eigenvalues of the matrix C on the diagonal. But I\ — I, hence the eigenvalues 
of Bj are the same as the eigenvalues of T (although with repetition). In other words, eigenvalues 
of T are eigenvalues of 3 as well. 

E.2 Eigenvalues of B 3 



Let u 



be a corresponding eigenvector of B3 (uj and u 2 are of size N x 1 each and as 

the eigenvalues of B3 are also eigenvalues of 3, we use Ag for an eigenvalue of B3). Hence, the 
standard eigenvalue relation B3U = Agfl requires the following equations to be satisfied: 



(a EE I + ps(3iA)ui + (ai E I + ps(3 2 A)u 2 = AgGi 

oceiUi + a n u 2 = \3VL2 



(21) 
(22) 
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Using Equation [22] we can compute ui in terms of u 2 as: 



ui 



a E i 



u 2 



(23) 



Substituting it back into Equation [21] we get: 

'Ag - OCn 



(cxeeI + pSPiA) 



ai E I + PsP2A u 2 = A 



Ag — a n 



which finally gives, 



a IE a E iI+ (ps(3i(A 3 - a n ) +Psp2a E i)A)u 2 



a E i 

Ag(Ag - a n )u 2 



u 2 



Au 2 



(a n + ocee)A 3 + aiia E E - ai E a E i 
PsPi(Ag - an) +PgP2a E i 



u 2 



(24) 



Again, Equation 24 tells us that either u 2 = or it is an eigenvector for A. But u 2 = => Ui = 
6 u — which is not possible. Thus Equation 24 is an eigenvalue equation for the adjacency 
matrix A and we are looking for solutions Ag and u 2 such that they satisfy it. Hence, 



(a n + a E E)Ag + ana E E - ai E a E i 



PsPi(Ag - a n ) +Psp2a E i 
where Aa is an eigenvalue of A. This finally gives 

Ag - Ag(a EE + a n + Ps0iA A ) + (a n a EE - a IE a E i + plA A (|3 1 a II - (3 2 a E i)) = (25) 

Thus we have a different quadratic equation (Q.£.) for each eigenvalue Aa of A. Each Q.E. gives 
us two eigenvalues (possibly repeated) of 3- 

So, finally, we can conclude the following lemma: 



Lemma 1 (Eigenvalues of 3)- Eigenvalues of 3 are given by the eigenvalues o/T (Equations 19 and 20 1 
and the roots of the Q.Es given by Equation^25\for each eigenvalue Aa of A. 



F Stability 

We require that all the eigenvalues of 3 to be less than 1 in absolute value (according to Theo- 
rem]^. From Lemma [T] we have two cases to handle in enforcing this: 

(CI) All the eigenvalues of T should be less than 1 in absolute value 

(C2) All the roots of the Q.Es given by Equation [25] for each eigenvalue A A of A should be less 
than 1 in absolute value 



F.l Case CI 

Note that this case depends only on the model as the matrix T is independent of the adjacency 
matrix A. But T is a stochastic matrix i.e. all the column sums are equal to 1 - consequently all its 
eigenvalues are less than 1 in absolute value. 



Lemma 2 (Stability CI). All eigenvalues of the matrix T (given by Equation 20* are less than 1 in 
absolute value. 
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F.2 Case C2 



As CI is always true, we need to only ensure case C2. We can prove here the following: 



Lemma 3 (Stability C2). All the roots of the Q.Es given by Equation 25 for each eigenvalue Aa of A are 
less than 1 in absolute value if. 

,( Pi(l-a„) + P2«Ei \ . 

AlPs 7^ > < 1 

V (1 - aii)[l - a EE ) - aiEati/ 



Proof Let x\ and r 2 be the roots of Equation 25 (y\ and r 2 can be real or complex depending on 
Aa). Then we want 

|ri| < 1 and |r 2 < 1 

T\ and r 2 are real As the roots are real, Aa is such that the discriminant D of the quadratic 
equation is greater than zero. In this situation: 

ItxI < 1 and |r 2 | < 1 =>■ r a e (-1,1) and r 2 e (-1,1) (26) 

From the theory of quadratic equations, it is well known (see e.g. \ V7\) that for real roots x\ 
and x 2 of a Q.E. f(x) = ax 2 + bx + c (with a > 0) to lie in the interval (—1,1) the following 
conditions must be true: 

a — c > 0, 
a - b + c > 0, 
a + b + c >0. 

Intuitively, the first condition forces the product of the roots to be less than 1 while the last two 
conditions state that value of f (x) at —1 and 1 should not be "too small". In our case, these then 
translate into: 

a n a EE - a IE a E i +p|A A Oia n - (3 2 a EI ) < 1 (27a) 
1 + oc EE + a n +ps(3iA A + ana EE - ai E a E i + ?* s \ A {&ictn - |3 2 a E i) > (27b) 
1 - a EE - a n -ps|3iA A + kijOee - aiEa E i +PsA A ((3ia n - (3 2 a E i) > (27c) 



Equations 27b and 27c can be written as: 



AaPs L + ) < 1 (28a) 

AKS V(l + a„)(l + a EE )-ai E a E iy 

AKS V(l-an)(l-a EE )-ai E a Ei ; 

respectively. The above equations should be true for any eigenvalue A a of A which makes T> > 0. 
Recall that we are considering only undirected graphs, hence A is a symmetric binary (0/1) 
square irreducible matrix. As a result firstly, all its eigenvalues are real. Secondly, from the 
Perron-Frobenius theorem |16J the algebraically largest eigenvalue Ai of A is a positive real num- 
ber and also has the largest magnitude among all eigenvalues. Hence if the above equations are 
true for Aa = Aj we are done. Now note that 

(1 + a n )(l + a EE ) - ai E a E i > (1 - a n )(l - oc EE ) - ai E a E i 

and that 

(3i(l - an) + (3 2 a E i > -|3i(l + an) + |3 2 a E i 
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In addition the L.H.S in both equations is positive under Aa = Aj. So Equation [28a] is always true 
if Equation 28b holds (under Aa = Ai) i.e. 



Aipj 



|3i(l - on) + P2QCE1 
;i - an)(l - a EE ) - ai E a E i 



< 1 



(29) 



As Ai is the largest eigenvalue both algebraically and in magnitude, under Equation 29 

(Plan - |3 2 a EI ) 



1 - a n oc EE + ai E a E i -psA A (Pioc n - |3 2 a E i) 

'(1 - an)(l - a EE ) - a iE a E i 



> l-a n a Et 



ai E a E i - 



(3i(l - a n ) + (3 2 a E i 
((1 - an) 2 + ai E a E i) |3i + a E i(2 - a n - a EE )|3 2 



(1 - a n )(3i + a E iP 2 



> 



.'. Equation |27a| is also true if Equation |29]holds. Thus the condition for the roots to be in (—1,1) 
when they are real is given simply by Equation '. 
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Ti and r 2 are complex In this case Aa is such that D < 0. Also as Equation 25 has real co- 
efficients, ri and r 2 are complex conjugate of each other and so |ri| = |r 2 | = V / Tv r r 2 . But the 
product of roots %i and x 2 of the equation ax 2 + bx + c = is equal to c/q. Hence we want to 
enforce c/a < 1. In our case it is 

anc* EE - ai E a E i + psA A (|3ian - (3 2 a E i) < 1 



which is exactly Equation 27a From the above analysis, we already know that it is true if Equa- 
tion 29 holds. So, for any eigenvalue Aa for which D < 0, the roots have magnitude less than 1 
given Equation 29 is true. 



Thus in both cases, whether roots are real or complex, Equation 29 is a sufficient condition for 
the roots to have magnitude less than 1. □ 

To re-cap we state Theorem [T] and then give its proof: 

Theorem (Super-model theorem - sufficient condition for stability). For virus propagation models 
which statisfy our general initial assumptions and for any arbitrary undirected graph with adjacency 
matrix A and largest eigenvalue Ai, the sufficient condition for stability is given by: 



where, s (the effective strength) is: 



s < 1 



Ai • C 



and C is a constant dependent on the model (given by Equation 30 K Hence, the tipping point is reached 
when s = 1. 

Proof Lemma|2]and Lemma[3]ensure cases CI and C2 and hence together with Lemma[T]impry 
that the eigenvalues of the Jacobian 3 of our general NLDS computed at the fixed point x are less 
than 1 in magnitude if Equation [29] is true. 

.'. using Theorem[2] our general NLDS is stable at its fixed point x if Equation[29]holds. Recall 
that x is the point when there no infected nodes in the system (Appendix [C| and that this is the 
fixed point whose stability conditions determine the epidemic threshold (Section [4||. 

.•. finally we can conclude Theorem [T] with 



C 



VPM 



Ps 



(3i(l - an) + |3 2 a E i 
(1 - an)(l - a EE ) - ai E a E i 



(30) 
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and the effective strength s = \ ■ Cvpm- The parameter Cvpm is a constant for a given propagation 
model while the only parameter involved from the underlying contact-network is Ai, the first 
eigenvalue of the adjacency matrix. □ 
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